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Abstract 



The evolution of a small-amplitude localized vortex disturbance in an unbounded 
shear flow with the linear velocity profile is investigated. Based on the exact solution 
of the initial problem, a revision is made of the theoretical approach (suggested by 
Levinski 199f and subsequently further developed in a series of other publications) 
in which the vortex evolution is described in terms of Fluid Impulse of the vortex 
"core" . Although the theoretical predictions obtained on the basis of that approach 
were excellently confirmed in subsequent experimental studies, its inconsistency is 
demonstrated in this study. 

According to this solution, the localized vortex increases slowly (as power-law 
with the time) and attains an almost "horizontal" orientation, unlike the previous 
theory (Levinski 1991) that predicts the more rapid growth and vortex orientation 
at the angle of 45° to the flow direction. On the other hand, just the rapid increase 
and the angle of 45° to the outer flow direction are characteristic for hairpin vortices 
observed in turbulent boundary layers or artificially synthesized vortices in laminar 
boundary layers. 

Thus the issue of adequate theoretical interpretation of the evolution of localized 
vortices is again on the agenda. The remaining part of the paper presents the first 
steps in the solving this problem. In particular, the dynamics of the total enstrophy 
of vortex as the measure of vortex intensity is followed. The dependence of vortex 
amplification on its initial orientation is investigated. On this base the validity of the 
old idea of Theodorsen (1952) on the predominant formation of the 45° vortices is 
discussed. Also the tensor of enstrophy distribition (TED) is defined and it is shown 
that it may serve an effective tool in describing of the vortex geometry. 

The linear stage of Gaussian vortex evolution presented here provides a very 
suitable base for testing of further numerical simulation of the nonlinear stage. 
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I. INTRODUCTION 



Two main types of coherent vortex structures, which were first identified using flow 
visualization technique in experiments reported by Kline et al. (1967), form the basis of 
the present views of the structure of a turbulent boundary layer. Thus the presence in 
the wall-bounded flow of streaks, along which the streamwise velocity is lower than the 
average velocity at the same distance from the wall, results from the rise of low-velocity 
fluid from near-wall layers induced by long-lived vortex structures. These latter represent 
a pair of counter-rotating vortices extended along the flow direction (Bakewell & Lumley 
1967; Smith & Schvartz 1983). 

Another phenomenon, widely observed in turbulent boundary layers and referred by 
Kline et al. (1967) as "bursting", results from the rapid evolution of the localized vortex 
having the shape of a hairpin. These vortices were found to be inclined at 45° to the flow 
direction (Head & Bandyopadhyay 1981), and their typical lifetime are about 5% of the 
streak lifetime. 

The evolution mechanisms of these well-organized vortex structures and their interaction 
have been subjects for study by an ever increasing number of researchers (see reviews by 
Robinson 1991 and Smith & Walker 1983). 

Despite the fact that both types of coherent vortices have an identical structure of the 
type of vortex dipole, their properties, and also formation and development mechanisms 
differ greatly. The slow evolution of near-wall vortices is reasonably well explained by 
the mechanism of algebraic growth suggested by Benney & Gustavsson (1981) and sub- 
sequently further developed in terms of the concept of optimal disturbances by Butler & 
Farrell (1992), Reddy & Henningson (1993) and Reshotko & Tumin (2001). On the other 
hand, the derivation of an adequate theoretical model describing the evolution of hairpin 
vortices is complicated by the high degree of vorticity localization in the core of the hairpin 
vortex and, hence, by the strong nonlinear character of their development from the outset. 
This is confirmed by a number of experiments where hairpin vortices were observed only at 
the transition of the threshold value of a certain parameter corresponding to the mechanism 
of their generation used in the experiment. Thus in experiments of Asai & Nishioka (1995) 
where hairpin vortices were generated by using acoustic disturbances, these vortex structures 
were observed only when the amplitude of the applied disturbance stood out above a certain 
critical value. In experiments of Malkiel, Levinski & Cohen (1999) the initial disturbance 
was created by employing suction of fluid through the holes in the wall. Here, as in the 
previous case, hairpin vortices were observed only at fluid suction rates exceeding a certain 
critical value. 

Based on the aforementioned factors, it is of interest to analyze the theoretical model 
describing the evolution of a nonlinear localized vortex disturbance in the external plane 
shear flow first suggested by Levinski (1991) and subsequently generalized to rotating flows 
(Levinski & Cohen 1995), flows of weakly conducting fluid in a magnetic field (Levinski, 
Rapoport & Cohen 1997) and to stratified flows (Levinski 2000). In this model the vorticity 
distribution is characterized by its fluid impulse integral defined as 




(1.1) 
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where r is the position vector, <jj(t, r) is the instantaneous field of vorticity disturbance, 
and the integration is done over the entire volume of fluid. Accordingly, the fluid impulse 
dynamics is described by the equation 



where the evolution equation for disturbed vorticity in a steady-state external velocity field 
(U) is obtained by applying the curl operator to the Navier-Stokes equation and by a 
subsequent substraction of the equation for undisturbed flow: 

— + (U ■ V) u> - (w • V) U - (ft ■ V) u + (u ■ V) u - (w • V) u = v Au>. (1.3a) 
at 

Here O = curl U, and u designates the disturbance-induced velocity field, 



Because the fluid impulse is invariant with respect to the self-induced motion of the vortex 
disturbance (Batchelor 1967), the nonlinear terms (underlined in equation (1.3a)) make a 
zero contribution to (1.2). This makes it possible to "linearize" the problem on evolution of 
a strongly nonlinear localized disturbance (for a more detailed description see the works by 
Levinski 1991 and Levinski & Cohen 1995, hereinafter L&LC). The fluid impulse integral 
has also an additional important property. Just by its definition (1.2), the fluid impulse 
describes both the increase in amplitude and the geometrical growth of the vortex. This is 
a highly important factor because the experimentally observed growth of hairpin vortices is 
not necessarily associated with the increase in vorticity amplitude. To cover such a scenario 
of a localized vortex evolution in terms of classical linear stability theory requires a very 
extensive analysis of amplitude changes for a great number of modes. 

The most important result, obtained on the basis of the fluid impulse approach, is the 
prediction of an exponential instability of a localized vortex disturbance in plane Couette flow 
(L&LC). This finding makes it possible to explain the experimentally observed formation 
and fast development of hairpin vortices in the boundary layer as resulting from a plane shear 
flow instability to localized vortices originating on the wall inhomogeneities. By applying 
this approach to circular Cuette flow (Malkiel, Levinski & Cohen 1999; Levinski & Cohen 
1995), it was possible to predict the growth of hairpin vortices within the range of basic flow 
parameters where the flow is known to be linearly stable. The criterion obtained for hairpin 
vortices growth was supported by the results of experiments performed by Malkiel, Levinski 
& Cohen (1999). 

On the other hand, as the procedure of deriving the closed evolution equation for fluid 
impulse, suggested in L&LC does not contain any formal limitations on the initial distur- 
bance amplitude, this result is in conflict with the results of classical linear stability theory 
for plane and circular Cuette flows (Drazin & Reid 1981, Dikii 1976). 

The objective of this paper is to analyze the evolution of a localized vortex disturbance 
in terms of linear stability theory on the basis of constructing a complete vorticity field. 
This approach is free from the deficiencies of the description of vortex evolution using the 
fluid impulse; unfortunately, however, it does not permit us to advance into the region of 





uj = curl it. 
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strong (nonlinear) vortices by analytical methods. On the other hand, the approach that is 
developed in this paper, is useful for analyzing the validity of the assumptions made in L&LC 
in the course of deriving the closed evolution equation for fluid impulse of the disturbance. 
Moreover, a knowledge of a complete vorticity field in the physical space for an arbitrary 
instant of time provides the basis for constructing the other integral characteristics of the 
disturbance vorticity field which can be useful in the analysis of the results of numerical 
simulations. 

This paper is organized as follows. 

In § 2 we subject to a critical analysis the fluid impulse concept suggested in L&LC, and 
on the basis of the exact solution of a linearized problem (for the time being, as a Fourier- 
representation) we demonstrate its invalidity for the present formulation of the problem (i.e. 
the problem of localized vortex evolution in the external flow with a lin ear velocity profile). 

In § 3, on the basis of the exact solution obtained above (in § 2), inverse Fourier-transform 
is used to construct the vorticity field in the physical space (within the linear approximation, 
of course). By considering an example where a so-called "Gaussian vortex" serves as the 
initial disturbance, the linear evolution of the vortex is studied for some particular cases of 
its orientation. It will be shown, in particular, that the symmetry properties of the basic 
equations forbid the formation of "hairpins" within the framework of a linear problem. 

In §4, we introduce the notion of total enstrophy of the vortex as the measure of its 
intensity. On this basis, we investigate the character of enhancement (attenuation) of the 
vortex depending on its initial orientation. 

In §5, we turn from the description of the vortex development based on its complete 
vorticity field to the description using a new integral characteristic, namely, the Tensor of 
Enstrophy Distribution (TED) which we introduce. This integral characteristic makes it 
possible to describe the vortex using only six independent parameters (and in the case of 
vortices symmetrical about the plane z — 0, even four parameters only). 

In § 6 we discuss the results obtained and the possible further directions of research. 

II. THE EVOLUTION OF THE FLUID IMPULSE. THE EXACT SOLUTION OF 
THE LINEAR INITIAL PROBLEM FOR THE VORTICITY FIELD IN THE 

K-SPACE 



It follows directly from the definition of fluid impulse integral (1.1) that it exists and is 
absolutely convergent only if 



Initially, a well-localized disturbance induces a velocity field possessing the asymptotic 
behavior (Batchelor 1967) 



A. Definition of the modified fluid impulse of vorticity 




where e > 0. 



(2.1) 




(2.2) 
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From substitution of (2.2) into (1.3a) it follows that the vortex field generated at an arbitrary 
time t > 0, has the asymptotic representation 

which does not satisfy the condition (2.1). 

To overcome this problem the approach based on vorticity separation procedure has 
been suggested by Levinski (1991). Accordingly, the vorticity field of the disturbance is 
subdivided on a closed vorticity field bounded the region directly adjacent to the initial 
vortex disturbance (u/), and vorticity field (w 11 ) that includes vortex tails generated in the 
process of the vortex disturbance development. Furthermore, it is assumed that uj 1 describes 
the evolution of the hairpin vortex, whereas uj 11 represents a vortex cloud, which has no 
substantial impact on the evolution of the concentrated vorticity. 

Since the subsequent analysis is based on the exact solution of linearized equations of 
vorticity dynamics and does not involve any additional assumptions, it is appropriate to 
introduce the concept of fluid impulse without recourse to the vorticity subdivision proce- 
dure. This will permit us, in particular, to analyze the validity of the procedure sussested 
by L&LC. 

For this purpose it is convenient to define the Modified Fluid Impulse (MFI), as follows 

pit) = lim \ [ r x u(t,r)dV. 

r<R 

The definition (2.4) is valid both in the case of the initial disturbance with an infinitely 
small amplitude as well as in a strongly nonlinear case. The only limitation is the local 
character of the disturbance at the initial instant of time, which corresponds to the absolute 
convergence of the fluid impulse integral for the initial distribution of vorticity. 

Since in the subsequent discussion the solution to the equations for vorticity dynamics 
is constructed in the Fourier-space, we shall use a Fourier-transform of the definition of the 
MFI. Namely, defining the Fourier-transform of the vorticity field as 

w(t,fc) = ( 27r r 3 / dVw(f,r) exp(-ifcr) 
the MFI (2.4) be represented as 

p(t) = ii(27r) 3 lim(V fc xa;(t,fc)), 

where the angle brackets correspond to averaging over the angles in the k-space: 

(•••) = — / cos Pdp d0 (•••), (2.7) 

47T J-tt/2 JO 

and (5 and (j> are the spherical angles in the k-space (with the axis y as the vertical axis, and 
the plane (xz) corresponding to (5 — 0): 

ki — k cos P cos 0, k 2 = k sin /3, k 3 = k cos f3 sin 0. (2.8) 
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(2.3) 



(2.4) 



(2.5) 



(2.6) 



It can be shown that the MFI defined by (2.4) is sufficiently "good" at first glance. 

Firstly, it does satisfy the necessary requirement of invariance with respect to the position 
of the center of a sphere. 

It should be noted, however, that in the present case where the vorticity decreases toward 
the periphery only as ~ M~ 4 (and the fluid impulse in the usual sense does not exist), the 
velocity field u(r) at large distances is no longer a potential one and, in particular, it no 
longer may be represented in the usual "dipole" form 

u = i cml {^) +0 ^ 

even if the MFI p serves as p. 

This means in particular that the MFI does not reflect at all the dipole structure of the 
disturbance vorticity field. As will be shown below, the plane of the vortex core localization, 
which can be described by the enstrophy distribution, L = |o;(r)| 2 , is not necessarily per- 
pendicular to the MFI direction, as should be in the case of the usual vortex dipole. (It can 
be shown, however, that even in this case the velocity again decreases toward the periphery 
in inverse proportion to the distance cubed, u ~ \r\ ~ 3 , as in the case of the potential velocity 
field induced by well- localized vortex). 

Secondly, it can be shown that the integral (2.4) exists at any instant of time and its 
value for sufficiently large values of R is independent on the value of R provided the MFI is 
well-defined at the initial instant of time. 

Indeed, by taking the time derivative of the expression (2.4) and substituting (1.3) into 
the right-hand side, we obtain 

= -\ hm J rx [(E7V) w - (wV) U - (OV) u] dV, (2.9) 

r<R 

Note that in (2.9) there are no contributions from the nonlinear and viscous terms (cf. (1.3)). 
The volumetric integral of these terms can be transformed to the integral over an infinite 
surface. The latter is zero by virtue of the asymptotic behavior of the vorticity u(r,t) and, 
accordingly, of the velocity field u(r,t) induced by it when R — > oo. 
Finally, the equation (2.9) may be transformed into form 

lim j. (210) 



where 



T J_ ^El 

■h- 2R £ ijk d 



j>^XiXjX m uj k dS 



1 dUi r 1 dU j r 

+ l~5 £ ij^ — f xix j x m u m db + —e ijk - — f xix k x m uidb. 
AR ox k Js AR J ox m Js 

Here eij k is the alternating tensor and usual summation convention is applied. By virtue 
of the asymptotic vorticity behavior, the limit of J is finite. This, together with (2.10), 
proves that if the MFI exists at the initial instant of time, then it exists also at an arbitrary 
instant of time. 
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Thus the MFI which we have just introduced seems, at first glance, a worthy replacement 
of the "fluid impulse of the core" introduced in L&LC for describing the evolution of a 
localized vortex since it does not require the usage of the vortex field subdivision procedure 
which has not been adequately justified in L&LC. 

It will be shown below, however, that any modification of the fluid impulse is unsatis- 
factory with regards to the capability to describe the structure of localized vortex. 

The left side of equation (2.10) together with the first two terms in its r.h.s. represents the 
evolution equation obtained in L&LC for the fluid impulse components constructed on the 
basis of a closed vorticity field a/, whereas the last term describes the specific contribution 
from the vortex tails. 

Note that the basic equation in L&LC is a linear equation in spite of the fact that 
it was derived from the exact nonlinear system of equations (1.3) without recourse to the 
linearization procedure. This means that the theory suggested in L&LC claims, in fact, a 
possibility of describing not only weak but also strong (nonlinear) vortices. In other words, 
this theory is insensitive to the vortex disturbance amplitude. 

On the one hand, this makes it extremely attractive, which, as a matter of fact, gave 
impetus to conduct a number of elegant experiments on its basis, and, on the other, if it is 
true, its predictions must remain valid for weak vortices as well. For weak vortices, however, 
there is a possibility of drastically simplifying the problem by performing a preliminary 
linearization of the initial system. This permits us to write the exact solution for the vorticity 
field and, on its basis, to check the validity of the theory. In particular, it is such a possibility 
of verifying the theory suggested in L&LC has stimulated this investigation. 

Thus, the contribution of vortex tails into (2.10) means that it is impossible to construct 
the closed equation describing the fluid impulse evolution (and from which the conclusion was 
drawn in L&LC about exponential instability) without one or another of vorticity subdivision 
methods. 

For that reason, below we make an attempt to give an alternative description to the 
evolution of a localized vortex without recourse to the evolution equation for fluid impulse. 
Using the solution obtained for the vorticity components we will also be able to calculate the 
fluid impulse and check to what extent the solution describing the fluid impulse evolution 
obtained in L&LC is consistent with what follows from this solution. 

B. The evolution of a localized vortex disturbance in the plane Couette flow 

1. Solving the initial problem for a localized disturbance 

It will be assumed that the basic flow has a linear velocity profile, U = (— fly, 0, 0) (so 
that its vorticity is 17 = (0,0, fi)). Note that this choice is not a loss of generality as a 
consequence of the assumption about the local character of the disturbance. It corresponds 
to the case where the characteristic size of the disturbance is much less than the characteristic 
size of variation of the basic velocity field. 

The suitable mathematical method for investigation of disturabance evolution in such 
flow was proposed by Lord Kelvin more than centure ago (Kelvin 1887, see also the more 
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recent publications based on this method, such as Craik & Criminale 1986, Criminale & 
Drazin 1990, Farrel & Ioannou 1993 and other), and we also apply it here. 

In this case the inviscid linearized equation for vorticity dynamics can be represented in 
the component-wise notation in Cartesian coordinates as 



n duJl 
ox 


n du3 

dx 


= 0, 


du 2 

OX 


n du2 

dz 


= 0, 


y n a 

OX 


n du3 

dz 


= o, 



(2.11) 



dt 
doo 2 
~dt 

d(jJ3 

dt 

where the subscripts "1", "2" and "3" correspond to the x-, y- and ^-components of the vec- 
tors, respectively. Upon introducing the dimensionless time r = Qt and Fourier-transforming 
the system of equations (2.12), we obtain 



du>i 
dr 


dk 2 


-ihu 3 = 


o, 


du 2 
dr 


dk 2 


- i k 3 u 2 = 


0, > 


du 3 
dr 


dk 2 


- i k 3 u 3 = 


o, 



(2.12) 



where Ui(t,k) represent the Fourier-transform of the disturbed velocity field components 
ut(r,t). 

The solution of the above problem on the evolution of the perturbation in the form 
of a single plane wave has been presented earlier in the work by Farrel & Ioannou, 1993 
(hereinafter F&I). Although the solution presented here below, in fact, is the same, we 
describe it derivation briefly to emphasize here on the vorticity components (instead of 
velocity components, as it done in F&I). We use, as in F&I, instead of the set of independent 
variables r, k±, k 2 and k 3 , a new set of independent variables r, ki, q and k 3 , where q = 
k 2 — k\T. We will consider uoi as a function of r, k±, q, k 3 , that is, uo,i = u>i(r, ki, q, k 3 ). Thus 
the new variable q is simply the initial value of the time-varying component of the wave 
vector k 2 . k 2 (t) = q + kvr. 

The new set of variables will be referred to as the Lagrangian coordinates in the k-space. 
The transition to the Lagrangian coordinates makes it possible to transform the system of 
equations (2.12) to a system of ordinary differential equations 





fci 


dr 


k 2 (r) 


duo 2 


Ht) 


dr 


k 2 {r) 


dw 3 


h 



dr k 2 (r) 



[k 2 {T)ui - kiu 2 \ = 0, 
[k 2 {r) uj 1 - kiu 2 ] + loi = 0, 
[k 2 [r) oj 1 - k x u 2 ] = 0, 



(2.13) 



where k 2 (r) = k\ + \k 2 {r)} 2 + k 2 , and d/dr = (d/dr) q = (d/dr) k2 + ki(d/dk 2 ) T . In deriving 
(2.13), we expressed also the Fourier-components of the velocity field in terms of Fourier- 
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k x uj(k^\ 

components of the vorticity field u(k) = i — . As a result we obtain the solution for 

dynamics of vorticity o?(r, k) as: 



A' 2 



Wl (r, fc) = ^(0, Q) - ^ r W2 (0, Q) + ^ V(r, fc) T(r, k), 

hp 



cu 2 (r, fc) = lu 2 (0, Q)--2- V(t, k) T(r, k), 



u 3 (r, k) = cu 3 (0, Q) - -y- r co 2 (0, Q) + 



V(r,k)T(r,k). 



(2.14) 



Here p — yk 2 + k 3 , Q — (q 1 , q 2 , q 3 ) = (fci, A: 3 ) is the initial value of the wave vector which 
in (2.14) must be expressed in terms of k: Q = (k 1: k 2 — fciT, fc 3 ), 

V (t, fc) = (1/p 2 ) [fc 3 Wi(0, Q) - hcu 3 (0, Q)] , T(r, fc) = arctan(£; 2 /p) - arctan[( k 2 -krf/p]. 

Including of the viscosity leads to addidional viscous factor in the expressions for uf 

i(r; k) -> ^(r; fe) exp(-(z//fi) jT fc 2 (r')dr') 



(see also F&I). 

In what follows, we shall use, as the initial vortex disturbance, the Gaussian vortex 



w (r = 0, r) = VF x /x, F = (tt 1/2 5)- 3 exp (-r 2 /5 2 ). 



(2.15) 



For numerical simulation of all hydrodynamic quantities (using program packages for 3-D 
hydrodynamics) it is often also necessary to specify the initial velocity field. It is readily 
calculated even for the initial isotropic function F(r): 



u = F(r) 



r {fxr)i H(r) 



t 1 ~ 



3r (fir) 



H(r) 



f 

Jo 



F(x) x 2 dx 



Thus we have in the ^-representation for the vortex of (2.15) 



k x fi) exp (— \k 2 S 2 ). 



(2.16) 



(2.17) 



Note that for the vortex (2.15) at the initial instant of time the usual fluid impulse p is also 
well defined and equal to fi. 

For our further purposes it is also very convenient to use the spherical coordinates in 
the Q-space: qi(— k±) = Q cos f3 cos 0, q 2 = Qsm{3 , q 3 (= k 3 ) = Qcos/9 o sin0, and 
— \k < p < |7r, < 4> < 2tt. In these variables we write finally for the Gaussian vortex 



Ui(r, k) 



pexp(-ig 2 J D 2 )O(r;/5 o ,0), 



(2.18) 
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Ci — A*3 tan @o — /i 2 sin — r cos 2 (/ii sin — /i 3 cos 0) — tan /3 cos 
C 2 = /ii sin - fi 3 cos + Q, 

C3 — ^2 cos — /ii tan flo — t sin cos (/ii sin — /i 3 cos 0) — tan /3 sin 
C2 = (Po ~ P) tan [tan /3q (//i cos + /x 3 sin 0) — /i 2 ] . 



Here 



D = D(t;Po,4>) = 6 { 1 + J^ [H-^cos / 9 o sin/5 o cos0+ §r 2 cos 2 A, cos 2 0]} i/2 , (2.20) 

and the "Reynolds number of vortex" Re is defined as Re = VL5 2 /v (so that in the inviscid 
case D — 8). The spherical angles (5 (in the fe-space) and f3 (in the Q-space) are related by 
the following expression (5 — (5 (r; /3q, 0) = arctan (tan fi Q + r cos 0). 



1/2 



(2.19) 



5. Evolution of the Modified Fluid Impulse of a localized disturbance 

As has already been discussed in the §1, the objective of this investigation was, in 
particular, to calculate the dynamics of fluid impulse of localized vortex in the external 
shear flow without using the vorticity subdivision procedure employed in L&LC. Indeed, the 
general solution of the vorticity dynamics equations, presented in §2.2.1, makes it possible 
to calculate the MFI defined by (2.4) at an arbitrary instant of time. 

For illustrative purposes we avail ourselves of the Gaussian vortex model (2.15) intro- 
duced above. It should be noted that the theory by L&LC is insensitive not only to the 
vortex disturbance amplitude but also to its form. Therefore, the initial disturbance can be 
chosen rather arbitrarily. Choosing it in the form a Gaussian vortex (2.15) optimizes calcu- 
lations substantially, but from the other side it is sufficiently representative model (see also 
notation in §6). The initial vortex is shown in figure 1, portraying the enstrophy isosurface 
o; 2 (0, r) = const that represents the surface of a torus. 
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FIGURES 



FIG. 1. Surface of constant enstrophy |w(0, r)| 2 = const of the initial Gaussian vortex. 

Note that for this model the size of the vortex disturbance core is specified by the value 
of the parameter 8, and the vortex plane and vortex lines that represent concentric circles 
are normal to the direction of the initial fluid impulse /x. 

Substitution of the expressions (2.18)-(2.19) for vorticity components into the expression 
(2.6) for the MFI gives: 

ftM = E%Wfe (2.2i) 

fin = 1 + | (sin 2 0tan/? o (A) -/?)) = 1 + \h, 

fl 12 = § (sin tan (j3 - #>)) = ± J 2 , h ( 2 .22) 

n 2 i = |r, n 22 = i, n 33 = 2 - n n = i - |/i . 

where the angle brackets correspond to averaging over angles in the fc-space, and (3q = 
(3o (r; [3, 0) = arctan (tan j3 — r cos 0). In the expanded form we have 

Px(t) =/i 1 + i/i 1 / 1 (t) + i/i 2 / 2 (t), p 2 (t) = // 2 + |a*i fit, p 3 (0 = A*s - \^h{t). (2.23) 

It is easy to obtain the asymptotic expressions for I x and J 2 for small (t <C l/|fi|) and large 
(f > l/|n|) times: 

hit) = l(nt) 2 + o((nt) 3 ), i 2 (t) = nt + o((nt) 3 ), \n\t<&i, (2.24) 
5 

and 

i 1 (t) = \n\t-3 + o(i/\n\t), i 2 (t) = 3[\n(\n\t) -0.6} + 0(1/ \n\t), \n\t^>i. (2.25) 

For small t we have from (2.23) and (2.24): 

Pi « Ml + 1/^2(^0) P2~M2 + |Ml(^), p3 = H3- (2-26) 

As would be expected, the MFI dynamics in the case of small times is determined by the first 
two terms in the evolution equation (2.10), in full agreement with the theory by L&LC. This 
is because the initial vortex (2.15) is well localized; therefore, the contribution J associated 
with the presence (at the early stage of evolution still very weak) vorticity "tails" at the 
periphery of the vortex, u> ~ r~ 4 , is vanishingly small. 
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In the case of larger times, however, the situation changes radically. The vorticity "tails" 
now become quite important and begin to affect the fluid impulse dynamics. 

This leads, in particular, to the fact that, according to (2.25), at large t the fluid impulse 
increases as a power law rather than exponentially, as was the case in L&LC. 

The fact that both the vorticity itself (including, of course, its "tails" produced in the 
course of evolution) and the MFI increase not more rapidly than as a power law, makes 
it possible to prove rigorously that there is no way to subdivide the disturbance vorticity 
field into two different components as it was suggested in L&LC. Namely, it is impossible 
to separate a localized vortex core, the fluid impulse of which grows exponentially, from the 
complete vorticity field (see §6 for more details). 

But it is the assumption about the possibility of such a separation that formed the basis 
of the approach suggested in L&LC. 

This means that the approach by L&LC is incorrect, in spite of a number of predic- 
tions obtained on this base which show an excellent agreements with experimental findings. 
Consequently, the problem of constructing an adequate theory describing the dynamics of 
localized vortices in shear flows becomes of current importance again. 

Nevertheless, it is interesting to point out that the inclination angle St (of the MFI vector 
to the positive direction of the x-axis tends to 45° with the time, in exactly the same way 
as does the fluid impulse p 1 constructed from the "core vorticity" in L&LC. 

Indeed, when \Q\t 3> I from (2.38) using (2.40) we have 

p 1 (t) = -^ 1 + fa \Q\t+\n2 \H\Vt\t) - 0.6], p 2 {t)=H2 + 

P3(0 = §A*3-|a*3 Mt, 

so that for vortices symmetric about the plane z = (113 = 0) we obtain: tan SI/ = 
p2(t)/pi(t) — > 1. Figure 2 shows the evolution of the quantities h(r) and ^(t) as well 
as of the inclination angle ^ for the case /x = (1, 0, 0). 




(2.27) 




FIG. 2. The evolution of I\{t) and hir) and of the inclination angle ^ of the modified fluid 
impulse for the vertically oriented (fii = 1 , fi2 = 0) symmetric (^3 = 0) Gaussian vortex. 
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However, in spite of the relatively good agreement between the orientation of the vortex 
plane observed in experiments and its orientation following from the description of the vortex 
using the MFI, it does not mean at all that the MFI is a more acceptable characteristic for 
describing localized vortices which should replace the "unfortunate" description using the 
"fluid impulse of the core" p 1 suggested previously in L&LC. 

It will be shown below that it is also possible to suggest some other methods of modifying 
the fluid impulse which will also be free from difficulties associated with the convergence of 
the corresponding integral at large distances, as is the just considered MFI; however, they 
will lead to a totally different scenario for the vortex geometry evolution. 

Actually, this would mean that the fluid impulse is not an adequate characteristic at all 
for describing the evolution of localized vortices. At least for the statement of the problem of 
localized vortex disturbance development in the external flow with the linear velocity profile 
accepted here (as well as in previous publications L&LC, Levinski, Rapoport & Cohen 1995, 
Malkiel, Levinski & Cohen 1999 and Levinski 2000). 

3. The evolution of the "Lagrangian" Modified Fluid Impulse of a localized disturbance 

In this section we introduce the concept of the fluid impulse of a selected ("colored") 
group of fluid particles and investigate its evolution. 

Let the initial position of a fluid particle that resides at the time r at a point with the 
coordinates r = (xi, x 2 , £ 3 ) be designated as r : r (r = 0) = r , where the components of 
the vector r be s 1 , s 2 and s 3 : r = (si,s 2 ,s 3 ). Then x 1 = s 1 — s 2 r, x 2 = s 2 , x 3 = s 3 , 
We now select a group of particles which at the initial instant of time are enclosed within 
a sphere of radius R and mentally paint it: ro = \J s 2 + s 2 + s 2 < R- We shall keep track 
on this painted group and calculate its fluid impulse. At subsequent instants of time, when 
r^O, the painted sphere will transform to an ellipsoid (x 1 + x 2 t) 2 + x\ + x\ < R 2 . 

Thus the fluid impulse of the painted group of particles is 

Pi(r, R) = \ e ijk / XjU k (r; r) dV, 

Jell 

where the integral is taken over the volume of the ellipsoid. By letting further R — > oo we 
obtaine for the LMFI components: pi = Tlijfij. Omitting the explicit expressions for tl ik , we 
present here only their asymptotic expressions for r ^> 1: 

n u w-i r + 2, fl 12 w -|r + |lnr- 0.414, fl 21 w -± r 2 + | r- 1 In r, 

n 22 ~ -| r 2 + \ t In r + 0.026 r, fl 33 w -\ r In r + 0.64 r. 

It is evident that for large r 

Pi ~ -\r (iii + 1^2), P2 ~ rpi ps -\t 2 (/ii + n 2 ), P3 ~ -|rlnr// 3 . 

We can now easily calculate the inclination angle ^ of the vector p = (]3i,]3 2 , ]3 3 ) in the 
(:q/)-plane for different orientation angles of the initial fluid impulse p(0)(= fx): tan^(r) = 

P2(r)/Pi(r) = [n 2 i/ii + n 22 // 2 ]/[n n //i + fi 12 fi 2 ]. 
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Figure 3 shows the evolution of *(r) for 8 orientations of fi: a n = (n — 1)45°. 




4 8 12 16 20 

FIG. 3. The evolution of the inclination angle * of the vector p. Numbers at the curves 
correspond to the initial orientation. 

It is evident that at large times the two-dimensional vector p = (p 1 ,p 2 ) is directed 
vertically. The vortices, the directions of the initial fluid impulse of which lie in figure 3 in 
the shaded and unshaded areas are directed at large r downward and upward respectively. 

It is interesting to note that over a sufficiently long time interval 2 < r < 6 the inclination 
angle is about 45° or 225° (according to the initial inclination angle a). 

Summary. Thus it is evident that two different methods of modifying the fluid impulse, 
MFI and LMFI, lead to two totally different results as regards their orientation at asymp- 
totically large times: 45° for MFI, and 90° for LMFI. It will be recalled that we associated 
intuitively the orientation of the vortex plane with the orientation of the fluid impulse vector 
p by assuming that, as for the usual dipole structure (such as in magnetostatics if we mean 
the analogy: u — > H, uj — > j, p — > m, where H is the magnetic field, j is electric current 
density, and m — \ J(r x j) dV is the magnetic dipole moment), this plane must simply be 
normal to the fluid impulse direction (L&LC) just as the plane of a ringlet with current is 
perpendicular to the dipole magnetic moment. 

It now becomes clear, however, that the fluid impulse in this problem just cannot describe 
adequately the vorticity distribution. By choosing in a different manner the form of the 
domain of integration, we can obtain for the same vorticity distribution not only an arbitrary 
time dependence of its fluid impulse but also an arbitrary inclination of the vortex plane. 

For that reason, there inevitably arises the problem of calculating a complete vorticity 
field. It is a fairly complicated numerical problem which is being solved to date (preliminary 
results of these calculations are presented in Suponitsky et al. , 2003, 2004), however within 
the linear approximation, it is actually solved (for single plane wave) by F&I and also in 
§ 2.2.1. It will now suffice to perform an inverse Fourier-transform and calculate the vorticity 
field in the physical space. This is done in § 3. 
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III. CALCULATING THE COMPLETE VORTICITY FIELD IN PHYSICAL 

SPACE 



We have 
where 



u)i(r) = J tOi(k) exp (i kr) d 3 k, 



(2vr) : 



■pd(k), 0(k) = 6(A0;r)exp(-iQ 2 5 2 ), 



(3.1) 
(3.2) 



and the quantities Q are specified by the expressions (2.19). 

The most compact method for evaluating the integral of (3.1) lies in passing from in- 
tegrating in the fc-space, to integrating in the Q-space (i.e. in the space of initial wave 
numbers). Introducing the spherical coordinates r , 9 and ip of the point r : 



r 



= \J (x~i + x 2 r) 2 + x\ + x~l, cos^o = 



COS(/?o 



Xi + X 2 T 



X2 



(3.3) 



yj (X 1 + X 2 T) 2 + X~l 

and also the angle 9 between the vectors Q and r : 

cos 6 = cos #o sin f3 + sin 9 cos f3 cos(0 — (p ). 



(3.4) 



we obtain 



Ui{r- r) = 



tt/2 



cos 2 /9 d(3 



x 



7T 5 /2 

r o cos0 o \/3 rQCOS 2 o 



J ^C,(/3o,0;r)(— — )(- 2 D2 y _ v D2 



) exp(- 



r-g cos 2 6 



(3.5) 



)■ 



In spite of the fact that the expression (3.5) is a sufficiently compact one, it still is very 
difficult for analysis, as it includes double integrals. For that reason, we have to carry out 
the subsequent analysis numerically. 

Note that the applicability of linear theory is limited by the condition |oj»| max \Q\. 
This means, in particular, that this condition must also be satisfied for the initial vortex, 
i.e. 

I^|max=|w(r = 5/V / 2)| = W4-H = 0- I 54^«fi, or /U «6.49fi5 4 . 

In order to investigate the linear evolution of the vortex, we calculated numerically the 
vorticity field distribution for fixed instants of time r by formula (3.5). Results are presented 
in Fig. 4 in the form of 3-D isosurfaces of absolute value of vorticity (or, that is the same, 
of the enstrophy density L) for fixed instants of time r: L(r; r) = |w(r; r)\ 2 = const. 
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T=0 T=1 T=2 T=3 




-1 -0.5 0.5 1 

x/6 

FIG. 4. The linear evolution of the Gaussian vortex: (a), (6) - horizontal and (c) - vertical. 
Isosurfaces of absolute value of vorticity are shown w(r;r) = 0.7u; max (r) = const, where cj max (r) 
is maximum (over volume) absolute value of vorticity at time r. 

From here on we shall confine ourselves to the case of symmetric initial vortices, 113 = 0. 
Furthermore, as is easy to understand, the vortex remains symmetric about the plane z = 
over the course of all subsequent evolution as well. Specifically, for the enstrophy density 
we have L(r; x, y, z) = L(r; x, y, —z). 

It is apparent from figure 4 that with the passage of time, the initial torus, corresponding 
to the Gaussian vortex, starts to rotate and deform and eventually turns into two symmetric 
"sausages" extended along the flow. (Note that only in this figure it is assumed that — Q, = 
dU/dy > in order to achieve a more usual perception of the vortex plane orientation.) 

These "sausages" could, in principle, serve as a source material for the hairpin legs. It 
can be shown, however, that within the framework of linear theory the sausages cannot turn 
into a hairpin through the formation of a bridge (a so-called hairpin head) near only one of 
the ends of the pair of legs. It turns out that this is forbidden by the symmetry properties 
of the basic equations! 

Indeed, it follows from the linearized set of equations (1.3) that if u>(0;r)= — u>(0; — r) 
and u(0; r) = it(0; —r), i.e. if all components of the initial vorticity change their sign with 
r — > —r and, accordingly, all components of the initial velocity do not alter their values 
in the case of the substitution r — > —r, then this symmetry property remains during the 
vortex (linear) evolution: 

u)(t; r) = — u>(t; — r) and u(t; r) — u(t; — r). (3.6) 
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Note that the Gaussian initial vortex does possess the aforementioned symmetry properties 

(3.6) . 

Consequently, for the enstrophy density L(t; r) of the initial Gaussian vortex at an 
arbitrary instant of time t we have L(t; x, y, z) = L(t; —x, —y, —z). In the case of a symmetric 
vortex, /i 3 = 0, we have an additional ^-symmetry: L(t; —x, —y, —z) = L(t; —x, —y, z). 
Consequently, we obtain L(t; x, y, z) = L(t; —x, —y, z). This means that a current enstrophy 
distribution must be invariant with respect to the simultaneous replacement x — > —x, y — > 
— y even for the same z. Specifically, directly "between the legs", i.e. in the plane z = 0, we 
have 

L(t;x,y,0)=L(t;-x,-y,0). (3.7) 

Hence we cannot obtain the "hairpin" in the course of the Gaussian vortex evolution, 
since the distribution of enstrophy in hairpin vortex does not have the symmetry property 

(3.7) . 

However, the situation is changed drastically if we include in consideration the nonlinear 
terms in equation (1.3) a). It is easy to see that the nonlinear terms (underlined in (1.3a)) 
totally destroy the symmetry properties of the linearized version of equation (1.3 a) and, 
hence, the prohibition for the hairpin in linear theory is removed] 

Therefore, (numerical) investigation of the nonlinear evolution stage of a localized vortex 
is strongly needed. Preliminary results of numerical calculations with strong vortices confirm 
the occurrence of hairpins at a definite evolution stage of the vortex (Suponitsky et al. , 2003, 
2004). 

IV. TOTAL ENSTROPHY OF A LOCALIZED VORTEX AND ITS GROWTH 

A. Calculating the total enstrophy of a localized vortex in inviscid flow 

In order to be able to describe the enhancement or attenuation of the vortex over the 
course of the evolution, we introduce, as one of its integral characteristics, the total enstro- 
phy: 

JHI£ = Ju;\r)dV (4.1) 

The total enstrophy C can serve as the measure of vortex intensity. Note that, at first 
glance, it seems more natural to take, as the measure of vortex intensity, its total energy £. 
However, an attempt to introduce its reasonable definition, like £ = /[(U + u) 2 — U 2 ] dV , 
runs into the same difficulty into which we ran in our attempt to introduce the fluid impulse. 
Indeed, since / (U • u) dV is divergent (remember that u ~ r -3 for large r), such a definition 
of the vortex energy cannot be recognized as correct. It is also easy to see that a similar 
introduction of the total enstrophy C is free from such difficulties because / u> dV = 0. 

The total enstrophy C depends on the time r as well as on the the initial fluid impulse 
[i. If the viscosity v is included, then C depends also on the Reynolds number Re. It 
is clear that finite viscosity effects are highly important for the problem of the maximal 
vortex enhancement (cf. also with F&I), especially in the connection with the problem of 
hairpin formation in the course of nonlinear evolution of initial weak vortex (which is partly 
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described by Suponitsky et al. , 2003, 2004 and will be considered in details in the following 
publication). 

Thus we have C = C(t, //). We can express the integral in terms of Fourier- variables: 

C= (2tt) 3 j \u(k)\ 2 d 3 k. (4.2) 

and use the corresponding expressions for Fourier-components of vorticity. We also introduce 
the normalized total enstrophy: 

£(t, M )=£(t,m)/£ (|m|), (4.3) 

where £o(|a*|) = (2n 3 )~ 1 / 2 5~ 5 |/x| 2 is the total enstrophy when r = 0. In the linear problem 
the normalized enstrophy £ depends only on the direction of \x. Consequently, it can be put 
without loss of generality that = 1. 

Finally, for the normalized total enstrophy we obtain £(r, fi) = %(r) /li/ij. 

B. Enhancement (attenuation) of the vortex and the relation to the hydrodynamic 

stability problem 

When r>lwe obtain 



l l3 r 2 + O (r In r), (4.4) 

9n 2 + 32 5 

hi = ~ 0.41954, Z12 = - « 0.55556, 

288 9 

0^2 _ An 3^2 

/ 22 = ^ ~ 1.35629, l 33 = — « 0.92529, 

36 32 

For the normalized enstrophy, when r>lwe have 

£(r)^a(^)r 2 , a(/z) = Zn yU 2 + 2 1 12 /ii/i 2 + Z22 ^2 + ^33 f4- ( 4 - 5 ) 

By analyzing the coefficient a(fi), we can readily find the orientation of the vector fj, corre- 
sponding to those initial vortices which will become the most enhanced at large times. By 
fixing \fjt\ = 1, we find that a is maximum when \x = (cosaoo, sin a^, 0), where 



«oo = 77 + 77 arctanf 2 ^ 12 ) « 65.07° 

Z 2 V ill ~~ '22 



I 2 
l 12 



1.6146. 



and is a = a max = \ + I22) + \J (hi ~ I22) 2 + 4ZJ 

Note also that the angle ~ 155.07° with a m i n ~ 0.1613. corresponds to the orientation 
of the least enhanced (at large times) vortices, i.e. the normalized enstrophy of the vortices, 
the initial orientation angle of which is directed along this direction, will be an order of 
magnitude smaller than the maximum one. 

Further, it is assumed again that [13 = 0. In this case it will suffice to describe the 
orientation of the initial vortex by only one inclination angle of its fluid impulse a: fi\ = 
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cos a, fi2 = sin a, /i 3 = 0. In other words, the angle a is the angle between the positive 
direction of the and the direction of the initial fluid impulse /x. 

For an arbitrary time r the maximum (over the whole range of initial orientation angles 
a) enhancement of the enstrophy corresponds to the initial inclination angle a{r) 

(4.6) 

A numerical calculation shows that a opt (r) increases from a opt (0) = 45° to a opt (oo) = 65.07° 



7T 1 

aopt(^) = ^ + -arctan 



-12 r 



£u(t)-£oo(t). 



60 



45 £ 1 , , , 1 1 , , T 

012345678 

FIG. 5. The inclination angle a op t(T) of the fluid impulse vector for which a maximum en- 
hancement of the enstrophy is reached by a given time r. 

Figure 6 shows the normalized enstrophy as a function of inclination angle a for four 
values of r r = 0; 0.5; 1.0 and 5.0. 

20 r , ^, 



'5 




FIG. 6. The Normalized enstrophy (radial coordinate) as a function of inclination angle of the 
initial fluid impulse a (angular coordinate) for four values of r: r = 0; 0.5; 1.0 and 5.0. 
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The calculations show also that at large r almost all enstrophy is concentrated in hor- 
izontal directions of vortex motion L\ ~ £3 = 0(r 2 £ 2 )- This fact is intuitively consistent 
with results of numerical calculations of the isosurfaces of enstrophy density (see §3, as 
well as § 5) which show that the localization plane of a weak vortex becomes horizontal 
asymptotically. 

As it follows from (4.5) for all possible orientation angles of the initial vortex we obtain 
a power-law (£ ~ r 2 ) growth of the total enstrophy (cf. with law of energy growth of the 
optimal inviscid excitations in F&I). 

We may regard such an asymptotic increase in enstrophy as the manifestation of an 
instability - the initial vortex starts to be enhanced infinitely (in linear theory at least). 

There is nothing surprising in this fact of a power-law enhancement of enstrophy. It is 
merely a reflection of another well-known fact, namely, that Couette inviscid flow is unstable 
with respect to 3-D wave disturbances, and this instability is just a power-law one, or, as it 
is also called, an 11 algebraic" instability (Ellingsen & Palm 1975). Recall that Couette flow 
is stable with respect to 2-D wave disturbances (k 3 = 0) altogether (that is, there is no not 
only the exponential instability but also even the algebraic instability). 

Note, that the last factor, together with the well-known Squire (1933) theorem stating 
that if the flow under consideration is stable with respect to 2-D disturbances it is nec- 
essarily stable with respect to 3-D disturbances), has long led to a paradox laying in the 
inconsistency between the experimental fact of the existence of turbulence in Couette flow 
and the absence of any instability in theory (for more detail see an excellent review by Hen- 
ningson, Gustavsson & Breuer 1994). This paradox was just resolved by the discovery of 
the power-law instability (u>i(r, k) ~ uj%(t, k) ~ r, uj2(t, k) ~ 1) of 3-D disturbances caused 
by a so-called lift-up effect (Landahl 1975) which wasn't taken into account in the proof of 
Squire's theorem. 

A localized vortex can be presented as a wave packet composed of 3-D wave disturbances. 
It is based on this that we have obtained here a power-law increase of the total enstrophy 
of the vortex. 

Precisely the same statement applies for the evolution of a localized vortex on the back- 
ground of Taylor-Couette (circular) flow considered by Malkiel, Levinski & Cohen 1999. 
Here also can not be an exponential growth of the localized vortex in the range of param- 
eters where the hydrodynamic stability theory predicts a stability. Now we have obtained 
the exact solution which describes the development of the weak vortex in the circular flow. 
It would be presented in separate paper. 

It is this fact that reflects the main conflict of the theory L&LC with well-known facts 
of the classical theory of hydrodynamic stability mentioned in § 1. 

Note also that including the viscosity is all the more unable to lead to an exponential 
growth of the vortex. On the contrary, the presence of viscosity leads to the fact that the 
vortex at some stage of its evolution can cease to increase and subsequently begin to be 
dissipated (see also F&I). Moreover, a sufficiently large viscosity unavoidably lead to the 
finiteness of the lifetime of a weak vortex. If, however, the viscosity is not too large and 
the initial amplitude of the vortex is not too small, the nonlinearity can come into play still 
before the vortex begins to be decay. This issue requires the further (numerical, of course) 
investigation. 
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C. Enhancement of the enstrophy and Theodorsen's idea of the predominant 

formation of 45-degree vortices 

As early as five decades ago Theodorsen (1952) came up with the hypothesis explaining 
why the 45-degree direction of orientation of horseshoe vortices, that is, the vortices whose 
plane is inclined at 45° to the basic flow direction, dominates in experiments. 

Since Theodorsen's (1952) idea is most clearly presented not in his work itself but in a 
later publication by Head & Bandyopadhyay (1981), we shall follow the presentation of this 
paper by adapting it to the present case of a linearized problem and an external flow with 
dU/dy = — £1 = const. 

We avail ourselves of the basic equation of the theory (1.3a) omitting the nonlinear (un- 
derlined) terms. A scalar multiplication of this equation by uj gives the equation describing 
the enstrophy dynamics of a fluid particle: 

j t (\u 2 ) = -nu^ + n (^^j +"( wAw ). ( 4 -7) 

where d/dt = d/dt + (UV) is the Lagrangian time derivative. The right-hand side of (4.7), 
with the exception of the viscous term, is the linearized "stretching term" ujujjdul/dxj, 
where the superscript "t" corresponds to the total (non-linearized) value of the physical 
quantity, u>* = Ct + oj, u f = U + u. 

By analyzing the first term on the right-hand side, it is easy to see that it is maximal 
when the two-dimensional vector uj 2 ) (at its fixed absolute value) is directed at the angle 
of 135° to the positive direction of the axis x, which corresponds to the angle of 45° measured 
from the direction of the mean velocity U in the upper half-space y > 0. (We are reminded 
that dU/dy < corresponds to positive values of Q. 

On the basis of this undeniable fact, Theodorsen arrived to the conclusion that the 
concentration of the vortex (its enstrophy) would also be maximal in the plane oriented at 
the angle of 45° to the flow, and this does correspond to experimental findings. 

Over the course of the past five decades Theodorsen's hypothesis has been repeatedly 
subjected to criticism from different standpoints. Here we want to discuss only two aspects, 
based on the just obtained (on the basis of the exact solution) results on the evolution of 
total enstrophy, results of calculations of the 3-D vorticity field presented in § 3, as well as on 
the results of calculations of the vortex localization plane inclination which will be described 
later in § 5. 

The first aspect implies that the conclusion about the greatest buildup rate of 45- 
degree vortices was drawn, strictly speaking, from analyzing the structure of only one of 
the terms on the right-hand side of (4.7). Our analysis shows, however, that the second 
term, Q(utdu/dz), that represents the other part of the "stretching term" responsible for 
the distortion of the flow velocity field caused by the vortex is also important. Although the 
integral contribution of this term to the total enstrophy is exactly zero at the initial instant 
of time, it becomes substantially larger with the time and can compete with the integral 
contribution of the first term. 

Thus the obviously true statement about the role of the stretching term made by 
Theodorsen can be applied, strictly speaking, only to the initial instants of time, r < 1, 



21 



when the distortion of the flow velocity field still can be neglected. Which one of the initial 
vortices will turn out to be the most enhanced for a sufficiently large time, in the course of 
which the direction of velocity is significantly changed, now becomes quite unclear from the 
reasoning presented. 

And the second aspect, which is of course associated with the first one, implies that at 
this point, taking into consideration the change of the velocity field orientation (actually 
neglected in Theodorsen's discussion), it becomes totally unobvious that the initial vortex 
with the optimal 45°-orientation, which at the initial instant of time was enhanced faster 
than all the others, would not change the orientation of its plane in the course of evolution. 

From the results presented in §4.2 (remember that they refer to the inviscid case) that 
are most instructively illustrated in figures 5 and 6, it follows that at small r the 45-degree 
vortices are indeed the strongest. However, at larger times, quite different vortices turn out to 
be most strongly enhanced, i.e. those for which at the initial instant of time the localization 
plane was more strongly pressed against the flow direction). In the limit r ^> 1 the vortices, 
which initially were inclined at 90° — a opt (r = oo) = 25° rather than 90° — a opt (r = 0) = 45° 
turn out to be the strongest. 

It should be noted at this point that, if the relatively small differences between the 
predicted angles of maximum enhancement are not taken into consideration, Theodorsen's 
hypothesis is plausible enough. 

As can be shown that including of the viscosity is also in favor of this hypothesis. At 
a finite Reynolds number, Re = Q5 2 /u, the difference between the angle of maximum 
enhancement and the angle of 45°, as predicted by Theodorsen, becomes still smaller. This 
can be most easily understood from the expression for total enstrophy at small r obtained 
with the including of viscosity (the details of its derivation are omitted): 

C(t, a; Re) = 1 + \r [ sin(2a) - — ] + O (r 2 ) . (4.8) 

It follows from (4.8) that initially the 45-degree (or, what is the same, 225-degree) vortices 
are the strongest. And the vortices whose plane is inclined at 135° are, on the contrary, the 
weakest (they are even weaker than the initial vortex, i.e. L < 1). It turns out that in the 
case of a sufficiently large viscosity the angles that are only very close to 45° are enhanced, 
and when Re = Re cr = 20 the only direction, 45°, is enhanced altogether. Calculation of the 
total enstrophy with the including of viscosity, shows that vortices of other directions, close 
to 45°, will start to be enhanced with the passage of time, however, the angles of maximum 
enhancement remain close around 45° during all time up to the beginning of the dissipation 
of the vortex. 

For illustration we presented in Fig. 7 the contours of £(r; a) for viscous case with 
Re = 20 and Re = 40. The shaded regions of the plane (r, a) correspond to C(t; a) > 1, i.e. 
to the enhancement of the vortex, and the unshaded areas correspond to C(r; a) < 1, that 
is, to the attenuation of the vortex. We see that if Re = 20 the vortex begin to dissipate for 
all a (except a = 45° and 225°), although with the growth of time some orientations near 
45° also (very weakly) enhance. 
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(a) Re = 20 




FIG. 7. Contours of the normalized enstrophy C(t; a) = const on the plane r — a for Re = 20 
and Re = 40. Unshaded areas correspond to the attenuation of the vortex, C < 1. 



It is interesting to note that this critical Reynolds number is very close to the value 
-Re m i n = 2ir 2 fa 19.7 obtained in F&I for growth of energy of optimal checkerboard excitations 
with k\ = 0. 

The situation with the second point concerning the inclination angle of the localization 
plane of the vortex is much less favorable. If we keep track on the evolution of the inclination 
angle of the enstrophy localization plane (vortex plane), we find that at large times this angle 
tends to zero (i.e. the plane of the vortex tends to become horizontal), rather than to 45°. 
This statement is illustrated in figure 4, as well as by calculations of this angle performed 
on the basis of the tensor of enstrophy distribution (TED) which are presented in § 5. 



V. THE TENSOR OF ENSTROPHY DISTRIBUTION AND VORTEX 

GEOMETRY 

We can often avoid an unwieldy description of the vortex by specifying its total vector 
field if we are able to introduce some integral characteristic of the vortex which (although 
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not reflecting, of course, in full measure the entire vector structure of the vortex) will permit 
its main geometrical parameters to be described at least roughly. 

For this purpose we avail ourselves of the analogy with those approaches which are used 
in electrostatics in describing the distribution of electric charge. This is customarily done 
by using so-called multipole moments. 

We now introduce the notion of the Tensor of Enstrophy Distribution, TED, which is 
essentially a usual quadrupole moment of the enstrophy distribution. 

We assume to use the following definition of the tensor (see, for example, the book by 
Levich 1969): 

(5.1) 



Ta 



, r XiX 



1^3 ■ 
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Note, that if the displacement of center of vorticity distribution takes place (i.e 
(/ dVu 2 Xj) I (J dVuj 2 ) 7^ 0, as it may be in the case of strong vortices) the definition of TED 
must be generalized: = J dVuj 2 (r) (x{ — Xi)(xj — Xj). 

As any symmetric tensor, it can be transformed to the principal ', where it has a 

diagonal form: 
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(5.2) 



Here Aj stands for the eigenvalues of the matrix T^, i.e. the solution of a characteristic 
equation Det | \Tij— A5y || = 0. Then the direction of one of the principal axes that corresponds 
to the smallest of these three values of A, is the direction which should be identified with 
a normal to the vortex plane. And the vortex itself is extended along the direction which 
corresponds to the largest value of A (see Fig. 8). 



short principal axis 



long principal axis 




FIG. 8. Illustration to the explanation of TED. 



It is also possible to introduce the notion of the size of the vortex a« along the corre- 
sponding principal axes x\: 

FF>~. 

(5.3) 
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Let the largest axes be referred to as a, the smallest as b, and let the third axis be c. 
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For illustration of the meaning of TED let us consider the initial Gaussian vortex. We 



have in this case: = K (2/z 2 £jj — /J,i[J,j) 



Ktij, where K 



l/(4 v / 2vr 3 / 2 (5 3 ). It is easy to 



see (transforming TED to the principal axes) that the direction along fi does correspond 
to the smallest of the three values of A, i.e. /z 2 , and in the plane which is normal to this 
direction, the eigenvalues are identical and are 2/z 2 . The size ratio along /x and in the plane 
perpendicular to \x is 1 : \/2. 

Thus the TED characterizes rather substantively the distribution of enstrophy for which 
we have an instructive visual idea from figure 1. 

Further, we again restrict our discussion to the case of a symmetric (about the plane 
z — 0) vortex, /i 3 = 0. In this case it is easy to see that the axis x 3 (axis z) remains 
one of the principal axes over the course of the entire evolution, and the TED has a more 
straightforward form 

A C 

C B . (5.4) 
D 

In this case the tensor is transformed to the principal axes by a simple rotation of the plane 
(xi,X2) around the axis x 3 (— x' 3 ) by an angle VP 



tan 2* = 



2C 
A-B' 



(5.5) 



and has in these axes a diagonal form. Let the axis x[ coincide with the shortest principal 
axis. Then in the new axes we obtain 
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(5.6) 



where S — y (A — B) 2 + 4C 2 and the angle between the positive direction of the axis x and 
the direction of a normal to the plane of the vortex (i.e. the axis x'{= x[)) 



* = iarctan(-4^) + i 



7T 



l + s), s = sign (A- B). 



(5.7) 



Next, using the notion of the TED introduced above, we can employ it to calculate the 
geometrical characteristics of the vortex and compare them with results that follow from 
calculations of the 3-D vorticity field by exact formula (3.5). 

Results of calculations of the vortex parameters, obtained on the basis of the TED, for 
four initial directions of the vector /x, a — 0°, 45°, 90°, 135° (where fi = (cos a, sina, 0)) 
are shown in Fig. 9. 
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1 2 3 4 5 

FIG. 9. The temporal dependence of the vortex parameters, 
initial values of the inclination angle a: (a) - a = 0°, (b) - 45°, (c) 
inclination angle $ as a function of r for these four initial values. 



The "axes" 
- 90°, (d) - 



a, b and c for the 
135° and (e) - the 



It is evident that at sufficiently large r the normal to the plane of the vortex is nearly 
vertical, m tt /2 — 1/t. 

In order to estimate the effectiveness of TED with respect to description of the vortex 
geometry, the parameters of the vortex geometry obtained from TED calculations were 
compared with results following from calculations of the complete vorticity field by exact 
formula (3.5). 

An exellent agreement of the inclination anlges ^ following from the TED with the 
actual inclinaton angles of the planes of enstrophy localisation was shown by Suponitsky, 
Cohen & Bar-Yoseph (2003, 2004). 

Thus it can be stated that the TED is a rather convenient and reliable integral char- 
acteristic for the description of the vortex dynamics. Of course, it is unable to describe 
the vector structure of the vortex, yet it can be used to obtain a sufficiently great deal of 
information about the vortex. 

The results on the orientation of vortex plane, obtained in this paragraph, unlike the 
findings concerning the evolution of total enstrophy, described in § 4 are virtually insensitive 
to the presence of viscosity. (It follows from our calculations of TED in viscous case, not 
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presented here.) 



VI. DISCUSSION 

As has been pointed out in the § 1, the motivation for this investigation was the analysis 
of a theoretical model suggested by Levinski (1991) for explaining the evolution mechanism 
of localized vortices observed in turbulent boundary layers. A key point in this model implies 
separating the complete vorticity field into the concentric vorticity with vortex lines enclosed 
within the region immediately surrounding the initial vortex disturbance, and the vorticity 
field associated with vortex "tails" which are produced in the process of evolution of the 
initial vortex disturbance. It should be noted that the possibility of such a separation is not 
strictly substantiated mathematically but is accepted in Levinski (1991) and in subsequent 
publications (Levinski & Cohen 1995; Levinski, Rapoport & Cohen 1997; Malkiel, Levinski 
& Cohen 1999 and Levinski 2000) as a physically justified hypothesis. The criterion of 
correctness of this approach comes from the agreement between results of theoretical analysis 
and experimental results. In particular, the predictions obtained on the basis of the model 
suggested in L&LC for rotating Couette flow were confirmed experimentally by Malkiel, 
Levinski & Cohen (1999). 

In this paper the hypothesis about the possibility of separating the vorticity is verified by 
constructing the complete vorticity field at an arbitrary instant of time for a small amplitude 
localized disturbance. The problem of the evolution of a weak localized disturbance is 
analyzed on the basis of exact solution for the external constant shear flow To ease the 
subsequent analysis of the vorticity field, the initial vortex disturbance was represented by 
"Gaussian vortex" (2.15) that specifies a very simple localized vortex, having the structure 
of a vortex dipole. 

It has been shown that in accordance with the classical stability theory results (see, for 
example, Dikii 1976), the vorticity amplitude increases not faster than it does as a power- 
law. This result contradicts the exponential growth of the fluid impulse obtained in L&LC 
for the "core" of the vortex disturbance identified in a special way. This could be accounted 
for by the fact that the generation of a new vorticity in the process of evolution of the 
vortex disturbance can lead to a fast increase of the "mass" of the vortex "core" . It is this 
phenomenon that is observed in visualizing vortex structures in turbulent boundary layers. 
Specifically there is a rapid growth (in the sense of the geometrical growth) of hairpin vortices 
which represent localized vortex dipoles. 

In order to analyze this possibility, we introduce the notion of the modified fluid impulse 
(MFI) defined as an integral of the dipole moment of vorticity over the infinite spherical 
volume. This definition coincides formally with the definition of the fluid impulse used in 
L&LC, but, unlike the latter, it is defined for the complete vorticity field. In doing this, we, 
using only the property of the localized character of the disturbance and without imposing 
constraints on its amplitude, show that if the MFI exists at the initial instant of time, 
then it exists also at any subsequent instant of time and does not depend on the particular 
coordinate system chosen. 

An analysis of the MFI behavior over large times t 3> shows that for any initial 

vortex orientation the MFI increases not faster than linearly with the time. 
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This result enables us to verify in a direct manner the hypothesis proposed in L&LC 
concerning the possibility of separation of the enclosed concentrated vorticity (localized 
vortex "core" ) from the complete disturbed vorticity field. Indeed, let us assume that such 
a separation is possible, that is, uj = uj 1 + oj 11 , where uj 1 describes the vortex "core", and 
(jj n describes the vortex "cloud" that includes all "tails" of the complete vorticity field. For 
each vorticity field, one can determine, in accordance with the expression (2.4), its modified 
fluid impulse, so that p = p 1 + p 11 Furthermore, as in L&LC, p 1 is a true fluid impulse. 
Accordingly, equation (2.10) that describes the dynamics of MFI defined for the complete 
vorticity field, breaks down into two equations 

dg = _ 1I dU 1 _ 1I dU 1 

dt " * P * d Xj * Pj dx^ { ' 



dpi i ~ii dUj ! u dUj_ 
dt dxi 2Pj dxi ' R-^'oo 



Pr -p- iPY ^+ lim MR)- (6.2) 



u 3 

Note that, according to the separation condition, the "tails" of the complete vorticity field 
make a contribution to the dynamics of p n only. 

An exponential growth of p 1 follows from equation (6.1). The fact that the sum p 1 + p n 
increases not faster than as a power- law, implies that p n also increases exponentially fast. 
It should be noted here that the term J in equation (6.2), describing the contribution from 
the vorticity "tails" to the dynamics of p H grow also not faster than as a power-law. Thus 
the main contribution to the MFI dynamics for the field uj 11 is made by the region that 
immediately surrounds the vortex "core" and, hence, the assumption about the possibility of 
separating the vortex core is invalid. 

On the other hand, the fluid impulse, defined for the complete vorticity field, cannot be 
an adequate characteristic of the evolution of a localized vortex. The formal reason is the fact 
that the volumetric integral involved in the definition of the fluid impulse is not absolutely 
convergent, and its value depends on the form of the integration domain when its size is 
made tend to infinity. In this paper this is illustrated by a comparison of the asymptotic 
values of the fluid impulse for two cases where the region of integration represents a spherical 
volume, first in Euler coordinates, and then in Lagrangian coordinates. 

In summarizing all attempts to describe the evolution of a localized vortex in the external 
shear flow, it can be stated that using the moments of the vorticity field in this problem is 
unjustified. 

In order to be able to describe the enhancement or attenuation of the vortex and the 
variation of its orientation the course of the evolution, we have analyzed the evolution of the 
total enstrophy of the vortex (4.1) and of the tensor of enstrophy distribution (5.1). They 
permit the evolution of the vorticity amplitude and the main geometrical characteristic of 
the vortex to be described by means of only a few independent parameters. 

In particular, the effectiveness of the description of the vortex on the basis of the tensor 
of enstrophy distribution (TED) can be demonstrated by comparing visual pictures of en- 
strophy density isosurfaces |w(r)| 2 = const, constructed on the basis of the exact solution 
for the complete vorticity field, with what follows from the description based on TED for 
the inclination angles of the plane of the vortex. Thus the TED is a reliable alternative 
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(to the fluid impulse) integral characteristic that enables an instructive representation of its 
evolution, instead of an unwieldy description using the complete vorticity field. 

Thus, the calculations done in this paper show that the linear evolution results in a pair of 
rollers lying in the horizontal plane and aligned along the flow. It is interesting to note that 
this result is in excellent agreement with the finding reported in F&I where the evolution of 
the form of energetic isosurfaces was calculated for optimal checkerboard perturbation. In 
spite of the large difference of the initial perturbations considered in our paper and in F&I, 
the outcome of their evolution turned out to be strikingly alike (see figures 9 and 10 in F&I). 
This indicates the universal character of the mechanisms giving rise to ordered structures in 
the course of the evolution of 3-D perturbations in shear flows, as declared in F&I, and also 
lends support to the idea that the particular form of the initial vortex selected in our paper 
is not very important. 

However, these results that follow from the exact solution of the evolution problem for 
the small amplitude localized vortex, contradict the known experimental facts obtained by 
visualizing hairpin vortices developing in turbulent boundary layers (Head & Bandyopadhyay 
1981) or artificially synthesized in laminar boundary layers (Acalar & Smith 1987 a,b). 

Moreover, as shown in § 3, the symmetry properties of the basic equations in the linear 
case, in principle, do not allow the formation of hairpin vortices. For that reason, (numerical) 
investigation of the nonlinear stage of evolution of a localized vortex is of utmost current 
importance. Preliminary results of numerical simulations with strong vortices confirm the 
occurrence of hairpins at a certain stage of vortex evolution (Suponitsky et al, 2003, 2004). 
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